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ABSTRACT 

Satellite kinematics can be used to probe the masses of dark matter haloes of central 
galaxies. In order to measure the kinematics with sufficient signal-to-noise, one uses 
the satellite galaxies of a large number of central galaxies stacked according to similar 
properties (e.g., luminosity). However, in general the relation between the luminosity of 
a central galaxy and the mass of its host halo will have non-zero scatter. Consequently, 
this stacking results in combining the kinematics of satellite galaxies in haloes of 
different masses, which complicates the interpretation of the data. In this paper we 
present an analytical framework to model satellite kinematics, properly accounting 
for this scatter and for various selection effects. We show that in the presence of 
scatter in the halo mass-luminosity relation, the commonly used velocity dispersion 
of satellite galaxies can not be used to infer a unique halo mass-luminosity relation. 
In particular, we demonstrate that there is a degeneracy between the mean and the 
scatter of the halo mass-luminosity relation. We present a new technique that can 
break this degeneracy, and which involves measuring the velocity dispersions using 
two different weighting schemes: host-weighting (each central galaxy gets the same 
weight) and satellite-weighting (each central galaxy gets a weight proportional to its 
number of satellites) . The ratio between the velocity dispersions obtained using these 
two weighting schemes is a strong function of the scatter in the halo mass-luminosity 
relation, and can thus be used to infer a unique relation between light and mass from 
the kinematics of satellite galaxies. 

Key words: galaxies: haloes — galaxies: kinematics and dynamics — galaxies: fun- 
damental parameters — galaxies: structure — methods: statistical 



1 INTRODUCTION 

According to the current paradigm, the mass of a dark mat- 
ter halo is believed to strongly influence the process of galaxy 
formation. Hence, a reliable determination of the masses of 
dark matter haloes can provide important constraints on 
the physics of galaxy formation. Numerous methods are 
available to probe the masse s of dark matter h aloes, in- 
cluding rotation curves (e.g., iRubin et al.l 119821) , grav ita- 
tional lens ing, either strong (e.g.. i Gavazzi et al.l l2007h or 
weak (e.g. iMandelbaum et al.l |2006|). X-ray emission from 
the hot intra-cluster medium fe.g.. iRvkoff et aUl2008l '), and 
a combin ed analysis of clustering and abund ances of galax- 
ies (e.g., lYang. Mo. fc van den Boschll2003n . In this paper 
we examine another powerful method, which relies on mea- 
suring the kinematics of satellite galaxies in order to infer 
the mass of the dark matter halo in which they orbit. 

In general, a dark matter halo hosts one central galaxy. 
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located at rest at the center of the halo's potential well, 
and a population of satellite galaxies orbiting the halo. 
The line-of-sight (hereafter los) velocity dispersion of these 
satellite galaxies reflects the depth of the halo's potential 
well, and is thus a measure for the halo's mass. In the 
case of rich galaxy clusters, the number of satellite galax- 
ies can be sufficient to properly sa mple the los velocity 
distributi on of its dark matter halo IjCarlberg et al.l Il996l : 
ICarlberg. Yee fc EUingsorJ Il997l ). In less massive haloes, 
however, the number of (detectable) satellite galaxies is gen- 
erally too small to obtain a reliable measure of the los veloc- 
ity dispersion. However, under the assumption that central 
galaxies with similar properties (i.e., luminosity) are hosted 
by haloes of similar masses, one can stack central galax- 
ies and combine the kinematic information of their associ- 
ated satellites to improve the statistics. Pioneering efforts in 
this d ir ection were made by Erickson, Gottesman & Hunteij 
(ll98it.lZaritsky et al.l (ll993friZaritskv fc Whitel l|l994) and 
IZaritskv et al.l ( 19971 ) . Although these studies were typically 
limited to samples of less than 100 satellites, they neverthe- 
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less sufficed to demonstrate the existence of extended dark 
matter haloes around (spiral) galaxies. 

More recently, with the advent of large homogeneous 
galaxy redshift surveys s uch as the Sloan Digital Sky Survey 
fSDSS: lYorket al.ll2000l ) an d the Two degree F ield Galaxy 
Redshift Survey (2dFGRS; IColless et al1l200lh . it has be- 
come possible to apply t his method to much larger sam- 
ples of satellite galaxies ( McKay et al. 2002; Prada ct al. 
20031: iBrainerd fcgpecianl (20031: lvan"dcn Bosch ct al...2004 : 



Conrov et all l2007l : iBecker et al.1 l2007h . These studies all 
found the los velocity dispersion of satellite galaxies, o-gat , to 
increase with the luminosity of the host (central) galaxy, Lc- 
This is in agreement with the expectation that more mas- 
sive haloes host more luminou s centrals. In a recent study, 
iNorberg. Frenk. fc Cold l)2008l 'l have shown that the quanti- 
tative discrepancies between these previous studies mainly 
owe to differences in the criteria used to select central hosts 
and their satellites. This underscores the necessity for a care- 
ful treatment of selection effects in order to extract reliable 
meiss estimates from satellite kinematics. 

Except for Ivan den Bosch et al.l (|200J), all previous 
studies have been extremely conservative in their selection 
of hosts and satellites. Consequently, despite the fact that 
the redshift surveys used contain well in excess of 100,000 
galaxies, the final samples only contained about 2000 — 3000 
satellite galaxies. This severely limits the statistical accu- 
racy of the velocity dispersion measurements as well as the 
dynamic range in luminosity of the central galaxies for which 
halo masses can be inferred. The main motivation for using 
strict selection criteria is to select only 'isolated' systems, 
with satellites that can be treated as tracer particles (i.e., 
their mass does not cause significant perturbations in the 
gravitational potential of their host galaxy). Let P(M\Lc) 
denote the conditional probability distribution that a central 
galaxy of luminosity Lc resides in a halo of mass M. If the 
scatter in P{M\Lc) is sufficiently small, preferentially select- 
ing 'isolated' systems should yield an unbiased estimate of 
{M){Lc), which is the first moment of P{M\Lc). However, 
very little is known about the actual amount of scatter in 
P{M\Lc) and different semi-analytical models for galaxy for- 
mation make significantly different predictions (see discus- 
sion in Norberg et al. 2008). If appreciable, the scatter will 
severely complicate the interpretation of satellite kinemat- 
ics, and may even cause a systematic bias (van den Bosch et 
al. 2004; More et al. 2008). Furthermore, even if the scatter 
is small, in practice, satellites of central galaxies stacked in 
finite bins of luminosity are used to measure the kinematics. 
If the satellite sample is small, one has to resort to relatively 
large bins in order to have sufficient signal-to-noise. There- 
fore, even if the distribution P{M\Lc) is relatively narrow, 
this still implies mixing the kinematics of haloes spanning a 
relatively large range in halo masses. 

In this paper we demonstrate that whenever the scatter 
in P{M\Lc) is non-negligible, the (Tsat(ic) inferred from the 
data has to be interpreted with great care. In particular, 
we demonstrate that there is a degeneracy between the first 
and second moments of P{M\Lc), in that two distributions 
with different {AI){Lc) and different scatter can give rise to 
the same (Tsat(ic). Therefore, a unique {A/)(Lc) cannot be 
inferred from satellite kinematics without a prior knowledge 
of the second moment of P{M\Lc). However, not all hope 
is lost. In fact, we demonstrate that by using two different 



methods to measure crsat(ic), one can actually break this 
degeneracy and thus constrain both the mean and the scat- 
ter of P{M\Lc). In this paper we introduce the methodology, 
and present the analytical framework required to interpret 
the data, taking account of the selection criteria used to 
identify the central host galaxies and their satellites. In More 
et al. (2008; hereafter Paper II) we apply this method to the 
SDSS to infer both the mean and the scatter of P{M\Lc), 
which we show to be in good agreement with the results 
obtained from clustering and galaxy-galaxy lensing analy- 
ses. In addition, in Paper II we demonstrate that (i) the 
scatter in P{M\Lc) can not be neglected, especially not at 
the bright end, and (ii) the strict isolation criteria generally 
used to select centrals and satellites result in a systematic 
underestimate of the actual (Af)(Lc). 

This paper is organized as follows. In Section [21 we 
present two different schemes to measure the velocity disper- 
sion, the satellite- weighting scheme and the host- weighting 
scheme. In Section [3] we present a toy model which serves 
as a basis for understanding the dependence of velocity dis- 
persion estimates on the different parameters of interest. In 
Sections [4] and [5] we refine our toy model by including selec- 
tion effects and by using a realistic halo occupation distri- 
bution (HOD) model for the central galaxies. We use these 
more realistic models to investigate how changes in the halo 
occupation statistics of central galaxies affect the velocity 
dispersion of satellite galaxies, and we demonstrate how the 
combination of satellite-weighting and host-weighting can 
be used to infer both the mean and the scatter of the mass- 
luminosity relation. We summarize our findings in Section [S] 
Throughout this paper M denotes the halo mass in units of 



2 WEIGHTING SCHEMES 

In order to estimate dynamical halo masses from satellite 
kinematics one generally proceeds as follows. Using a sample 
of satellite galaxies, one determines the distribution P{AV), 
where AV is the difference in the line-of-sight velocity of a 
satellite galaxy and its corresponding central host galaxy. 
The scatter in the distribution P{AV) (hereafter the veloc- 
ity dispersion), is then considered to be an estimator of the 
depth of the potential well in which the satellites orbit, and 
hence of the halo mass associated with the central. In order 
to measure the velocity dispersion as a function of central 
galaxy luminosity, (Tsat(ic), with sufficient signal-to- noise, 
one has to combine the los velocity information of satellites 
which belong to centrals of the same luminosity, Lc. This 
procedure is infiuenced by two effects, namely mass-mixing 
and satellite-weighting, which we now discuss in turn. 

Mass- mixing refers to combining the kinematics of satel- 
lites within haloes of different masses. The mass-luminosity 
relation (hereafter MLR) of central galaxies can have an 
appreciable scatter, i.e., the conditional probability distri- 
bution P{AI\Lc) is not guaranteed to be narrow. In this 
case, the satellites used to measure crsat(ic) reside in halo 
masses drawn from this distribution, and (Tsat(ic) has to be 
interpreted as an average over P{M\Lc). 

In most studies to date, the technique used to measure 
o"sat(ic) implies satellite weighting. This can be elucidated 
as follows. Let us assume that one stacks A'^c central galaxies. 
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and that the j^^ central has Nj sateUites. The total number 
of satellites A^'sat is given by ^ .^-^ Nj. Let AVij denote the 
los velocity difference between the i"' satellite and its cen- 
tral galaxy j. The average velocity dispersion of the stacked 
system, ctsw, is such that 
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Here aj is the velocity dispersion in the halo of the j^^ cen- 
tral galaxy. The velocity dispersion measured in this way is 
clearly a satellite- weighted average of the velocity dispersion 
aj around each central galaxy q Although not necessarily 
directly using Eq. ([T}, most previous studies have adopted 
this satellite- weighting scheme (McKay et al. 2002; Brainerd 
& Specian 2003; Prada et al. 2003; Norberg et al. 2008). 

In principle, the satellite- weighting can be undone by in- 
troducing a weight Wij = 1/Nj f or each satellite-central pair 
in the los vel ocity distribution (Ivan den Bosch et al.ll2004l : 
IConrov et al.l l2007l. The resulting host-weighted average ve- 
locity dispersion, (Thw, is such that 
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and it gives each halo an equal weight. 

Consider a sample of central and satellite galaxies with 
luminosities L > Lmin- The velocity dispersions in the 
satellite-weighting and host-weighting schemes can be an- 
alytically expressed as follows: 



^sw{Lc) = 



(^hw{Lc) 



/p°° P(MiLc) {JVeat)M (apM dM 
/;°P(M|L,)(iV,at)MdM 

/p°°P(M|Lc){aL)MdM 
/g°° P(M|Lc) dM 



(3) 



(4) 



Here (Nsat) j^.j denotes the average number of satellites with 
L > Lmin in a halo of mass M, and {o's^t)M is the square 
of the los velocity dispersion of satellites averaged over the 
entire halo. 

Consider a MLR of central galaxies that has no scat- 
ter, i.e. P{M\Lc) = 5{M - Mo), where Mo is the halo 
mass for a galaxy with luminosity Lc. In this case both 
schemes give an equal measure of the velocity dispersion, i.e., 
ag„ = a^^ = {<y1ax)Mo- Most studies to date have assumed 
the scatter in P[M\Lc) to be negligible, and simply inferred 
an average MLR, Mo(I/c) using a^^(Lc) = {a^i,t}Ma (McKay 
et al. 2002; Brainerd & Specian 2003; P rada et al. 2003; Nor- 
berg e t al. 2008). However, as shown in lvan den Bosch et al.l 
(|2004r ). and as evident from the above equations ([3} and (|4]), 
whenever the scatter in P{M\Lc) is non-negligible, a^^{Lc) 
and a^^{Lc) can differ significantljo (see also Paper II). 

In this paper, we show that ignoring the scatter in the 
MLR of central galaxies can result in appreciable errors in 



^ Note that the velocity dispersion is always averaged in quadra- 
ture. 

^ Note that a^^ ^ ct^^ is a sufficient but not a necessary condi- 
tion to indicate the presence of scatter in P{M\Lc); after all, if 
(Afsat)jvj does not depend on mass then a^^ = a^ independent 
of the amount of scatter. 



the inferred mean relation between mass and luminosity. We 
show, though, that these problems can be avoided by simul- 
taneously modeling it^„{Lc) and a^„{Lc). In particular, we 
demonstrate that the ratio of these two quantities can be 
used to determine the actual scatter in the MLR of central 
galaxies. 



3 TOY MODEL 

In the previous section, we have shown that both a-s^{Lc) 
and (T^„(Lc) can be analytically expressed in terms of 
the probability function, P{AI\Lc), the satellite occupation, 
{Ns!it)]^,j, and the kinematics of the satellite galaxies within a 
halo of mass M specified by {(Ts&t)M- In fact, the inversion of 
equations ((3} and Q presents an opportunity to constrain 
P(A/|Lc) using the observable (Tg^ and o'^^- In this section 
we use a simple toy model to demonstrate that the combi- 
nation of (Ts,v and a^^^ can be used to constrain the first two 
moments (i.e., the mean and the scatter) of P(M|Lc). 

For convenience, let us assume that P(M|Lc) is a log- 
normal distribution 



P(M|Lc)dM: 
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Here Mo is a characteristic mass scale which obeys 
lnMo= P(A/|Lc)lnMdM = (InM), 

"'0 

and fffn jvj reflects the scatter in halo mass at a fixed central 
luminosity and is given by 

/•oo 

crfnAf= / P(Af|Lc)(lnM-lnMo)'dM. 



(6) 
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In addition, let us assume that both {o'sa.t)M and 
(Nsat) ]^,i are simple power laws, 
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with a and 13 two constants, TV the average number of satel- 
lites in a halo of mass lO'^^ /i'^Mq, and S the corresponding 
los velocity dispersion. 

Substituting Eqs. ^, (O and (|9j in Eqs. Q and Q 
yields 



'^sw(^c) = 5' (^ j exp 



(^h^Lc) = S [-^j exp 



^2 fl2 
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(10) 

(11) 



The velocity dispersions CTsw(ic) and crhw(ic) depend on 
both Mo and ainM, elucidating the degeneracy between the 
mean mass Mo(I/c) and the scatter o"inA/(ic) of the distri- 
bution P(M|Lc). In particular, if only asw(i'c) or o-hw(^c) is 
measured, one cannot deduce A4o{Lc) without having an in- 
dependent knowledge of the scatter crinM(ic). However, the 
latter can be inferred from the ratio of the satellite- weighted 
to the host- weighted velocity dispersion. In particular, in the 
case of our toy model. 
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O-ln A/ = ^ in 



(12) 



Thus, by measuring both as„{Lc) and (Thw(Jjc) one can 
determine both Mo{Lc) and its scatter (TinA/(ic), pro- 
vided that the constants a and /3 are known. Since viri- 
alized dark matter haloes all have the s ame average density 
within their virial radii, / 3 = 2/3 (e.g. iKlvpin et al.lll999l : 
Ivan den Bosch et al.ll2004l ). Previous studies have obtained 
constraints on a that cover the range 0.7 < , a <, 1.1 (e.g. 
Yang et al.l |2005|: Ivan den Bosch et al.l l2007l : IXinker et all 



20071 : 1 Yang et al 



20071 1 . Since ai^ji^ x a , this uncertainty 



directly translates into an uncertainty of the inferred scat- 
ter. Therefore, in Paper II, we do not use the constraints on 
a available in the literature to infer the mean and scatter of 
the MLR from real data. Instead, we treat a as a free pa- 
rameter and use the average number of observed satellites 
as a function of the luminosity of central as an additional 
constraint. 



4 SELECTION EFFECTS 

The toy model presented in the previous section illus- 
trates that measurements of the satellite- weighted and host- 
weighted kinematics of satellite galaxies can be used to in- 
fer the mean and scatter of the MLR of central galaxies, 
P{M\Lc). However, in practice one first needs a method 
to select central galaxies and satellites from a galaxy red- 
shift survey. In general, central galaxies are selected to be 
the brightest galaxy in some cylindrical volume in redshift 
space, and satellite galaxies are defined as those galaxies 
that are fainter than the central by a certain amount and 
located within a cylindrical volume centered on the central. 
In this section we show how these selection criteria impact 
on a^.^ and (t^^, and how this can be accounted for in the 
analysis. 

No selection criterion is perfect, and some galaxies will 
be selected as centrals, while in reality they are satellites 
(hereafter 'false centrals'). In addition, some galaxies will 
be selected as satellites of a certain central, while in reality 
they do not reside in the same halo as the central (here- 
after 'interlopers'). The selection criteria have to be tuned 
in order to minimize the impact of these false centrals and 
interlopers. Here we make the assumption that interlopers 
can be corrected for, and that the impact of false centrals 
is negligible. Using mock galaxy redshift surveys, van den 
Bosch et al. (2004) have shown that one can devise adap- 
tive, iterative selection criteria that justify these assump- 
tions (see also Paper II). Here we focus on the impact of 
these iterative selection criteria on the satellite kinematics 
in the absence of interlopers and false centrals. Our analyt- 
ic al treatment for select i on effe cts follows the one presented 
in Ivan den Bosch et al.l (|200J) except for the inclusion of 
an extra selection effect. For completeness, we outline this 
treatment below. 

In general, satellite galaxies are selected to lie within a 
cylindrical volume centered on its central galaxy, and spec- 
ified by Rp < Rs and \AV\ < (AT/)s. Here Rp is the physi- 
cal separation from the central galaxy projected on the sky 
and AV is the los velocity difference between a satellite 
and its central. Usually, (A1/)s is chosen sufficiently large. 



so that it does not exclude true satellites from being se- 
lected . However, in the adaptive , iterative selection crite- 
ria of Ivan den Bosch et al.l l|2004 ). which we will use in our 
companion paper II, the aperture radius is tuned so that 
Rs — 0.375 rvir, where rvir is the virial radius of the dark 
matter halo hosting the central-satellite pair. This means 
that (A'aat)jvf and {as!^t)M in Eqs. (O and (|4)| need to be 
replaced by (A'^3at)ap,M and {crsat)ap,Af , respectively. Here 
{Nsa,t)ap,M is the average number of satellites in a halo of 
mass M that lie within the aperture, and {cr^at)ap,M is the 
square of the los velocity dispersion of satellite galaxies av- 
eraged over the aperture. 

The number of satellites present within the aperture, 
{A'sat)ap,Af , is related to the number of satellites given by 
the halo occupation statistics, (Nsat) j^j, via 



/cut (A''sat)j,j if -Rs < r-vir 
{Nsa.t)nj if -Rs ^ r-vir 



{iVsat)ap,A/ = <i f-^;'=-'« :;;;=-■- ^^g^ 



with 

Jcu 



Here nsa,t{r\M) is the number density distribution of satel- 
lites within a halo of mass Af , which is normalized so that 



(iVs, 



t)M = 47r / nsat(r|M)r^dr. 
'o 



(15) 



Under the assumption that the satellites are in virial 
equilibrium within the dark matter halo, and that the ve- 
locity dispersion of satellite galaxies within a given halo is 
isotropic, the los velocity dispersion of satellites within the 
cylindrical aperture of radius _Rs is given by 



\''"sat/ap,M 



47r 



{A'sat)ap,Af 



dRR 



nsat(r!Af)aL(riM) 



r dr 



V~r 



i?2 



(16) 



Here asat(f |Af) is the local, one-dimensional velocity disper- 
sion which is related to the potential ^ of the dark matter 
halo via the Jeans equation 



2 I 



\M) 



nsat(r|M) 
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risat(r'|M)--(r'|A/)dr'. (17) 
or' 



The radial derivative of the potential ^ represents the radial 
force and is given by 



d^ , ,,,, AttG 
or r2 



p{r \M) r dr , 



(18) 



with p[r\M) the density distribution of a dark matter halo of 
mass M. The assumptions of virial equilibrium and orbital 
isotropy are supported by results from numerical simulations 
which show that dark matter subhaloes (and hence satel- 
lite galaxies) are in a steady state equilibrium within the 
halo and that t heir orbits are nearly isotropic at l east in the 
centra l regions (|Diemand. Moore, fc Stadell 120041 ). Further- 
more, Ivan den Bosch et al.l (200J) have demonstrated that 
anisotropy has a negligible impact on the average velocity 
dispersion within the selection aperture. 

Finally, there is one other effect of the selection cri- 
teria to be accounted for w hich has not been considered in 
Ivan den Bosch et al.l (|200J). When selecting central-satellite 
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pairs, only those centrals are selected with at least one 
satellite inside the search aperture. This has an impact on 
the host-weighted velocity dispersions that needs to be ac- 
counted for. The probability that a halo of mass M, which 
on average hosts {Nsat)a.p,M satellites within the aperture 
Ra, has A'sat ^ 1 satellites within the aperture, is given by 



P(iVsat > 1) 



1 - P(iV,at = 0) 

1 - exp[-{7Vsat)ap,M] 

P({Afsat)ap,JVl). 



(19) 



Here, for the second equality, we have assumed Poisson 
statistics for the satellite occupation numbers. Note that, in 
the satellite- weighting scheme, haloes that have zero satel- 
lites, by definition, get zero weight. Therefore only the host- 
weighted velocity dispersions need to be corrected for this 
effect. 

Thus, in light of the selection effects, Eqs. (|3]) and Q 
become 



.{Lc 



and 



<^hw{Lc) = 



/p°° P(M|Lc) (JVsat)ap,M {aLt).p,M dM 
/g°°P(M|Lc){iVBat)ap,MdM 

/f P(M|Lc) P({7Vsat)ap,M) {aL)ap,M dM 
/g°° P(Af ILc) P({iVsat)ap,A/) dM 



(20) 



(21) 



Note that P({A'"sat)ap,M) ^ {A''sat)ap,Af whcU {iVsat)ap,M -^ 

0. This implies that |crsw — fhwl -^ for faint centrals (i.e. 
when I/c becomes comparable to the luminosity limit of the 
survey) . 



5 MORE REALISTIC MODELS 

Using the methodology described above, we now illustrate 
how satellite kinematics can be used to constrain the mean 
and the scatter of the MLR of central galaxies, P{AI\Lc). We 
improve upon the toy model described in Section [3] by con- 
sidering a realistic model for the halo occupation statistics 
and take the impact of selection criteria into account. 

As is evident from the discussion in the previous sec- 
tion, calculating (Ta„(Lc) and a^^{Lc) requires the following 
input: 

• the density distributions of dark matter haloes, p{r\M) 

• the number density distribution of satellites, nsa,t{r\M) 

• the halo occupation statistics of centrals, P(Af |Lc). 

We assume that dark mat ter haloes follow the NFW 
iJNavarro. Frenk. fc Whitelll997^ density distribution 



p{r\M) = 



M 



(^"■('-S 



(22) 



Here, r^ is a characteristic scale radius, c — r^^-frs is the 
halo's concentration parameter, and 



IJ.{x) = ln(l + x) — 



1+x 



(23) 



Throughout we use the relation between c and M given by 
iMaccidet ahl (|2007| ). 

We assume that satellite galaxies are spatially unbiased 
with respect to the dark matter particles, so that their num- 
ber density distribution is given by Eq. (|22|l with M replaced 
by (Nsa^t) J^J■ Note that this is a fairly simplistic assumption. 
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Figure 1. The satellite- weighted (ctsw) and host- weighted (crhw) 
velocity dispersions of satellite galaxies for model Gl. Note that 
o"sw(ic) > ChwC^c) at the bright end, indicating that the MLR of 
central galaxies, P(M|Lc), has anon-negligible amount of scatter. 



We address the issue of potential spatial antibias of satellite 
galaxies in Paper fl. 

Substituting p{r\M) and nsat(r-jAi") in Eqs. HH and 
(fTT)) gives 
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(^sat 



r-jA//) 



fi{c) 



(ao 



/>■= 



jj.{x) dx 
JJlTxj' 



(24) 



where Kir = (GAf/rvir)^ is the circular velocity at rvir. 

The final ingredient is a realistic model for the halo oc- 
cupation statistics of centrals and satellites. To that extent, 
we use the conditional luminosity function (CLE) presented 
in Cacciato et al. (2008). The CLE, denoted by ^{L\M)dL, 
specifies the average number of galaxies with luminosities in 
the range L±dP/2 that reside in a halo of mass M, and is ex- 
plicitly written as the sum of the contributions due to central 
and satellite galaxies, i.e. $(L|Af) = $c{L\M) + $s(P|A/). 
Erom this CLE, the probability distribution P(A/|Lc) fol- 
lows from Bayes' theorem according to 

$c(Pc|Af)n(Af) 



P(Af!Lc) = 



/J"$c(Lc|Af)n(Af)dAf ' 



(25) 



with n{M) the halo mass function, while the average number 
of satellites with L ^ Lmin in a halo of mass A/ is given by 



(^^sat)M = 



$s(P|Af)dL. 



(26) 



The parametric forms for ^c(L\AI) and $s(L\M) are 
motiv ated by the results of lYang. Mo fc van den BoschI 
(|2008l . hereafter YMB08), who determined th e CLE from 
the SDSS group catalogue of lYang et al.l l|2007l 'l. In particu- 
lar, "I>c(P|Af) is assumed to follow a log-normal distribution 
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Figure 2. Comparison of three models with different HODs for the central galaxies. In all panels the solid line corresponds to model Gl, 
the dotted line to model G2 and the dashed line to model G3 (see Table 1 for the parameters). Panels (a), (b) and (c) show (logLc)(A/), 
(IogM)(Lc) and CiQgjvf (^c)) respectively. Panels (d) and (e) show the predicted satellite-weighted and host-weighted velocity dispersions 
as function of luminosity, and panel (f ) shows the logarithm of the ratio between CTsw and crij^,, . See text for a detailed discussion. 



Table 1. Different models for the HOD of centrals 



Model 



""logL 



71 



72 



u 



Ml 



Gl 
G2 
G3 



0.14 
0.25 
0.14 



3.27 
3.27 
1.80 



0.25 
0.25 
0.40 



9.94 
9.94 
9.80 



11.07 
11.07 
11.46 



Three different models describing the MLR of centrals used to 
predict crsw(ic) and crijw(Lc). 



<E>c(i|M)di = 



loge 



2lT CTlog L 



exp 



log(L/L* 



v/2< 



Clog L 



¥,(... 



with (Tiog L a free parameter that we take to be independent 
of halo mass, and 



LUM) 



(M/Mi)^i 



[1 + (M/Mi)]^ 



(28) 



which has four additional free parameters: two slopes, 71 
and 72, a characteristic halo mass, Mi, and a normaliza- 
tion, Lo- Note that, L* ex M^i for M < Mi and L* oc M^^ 
for M ^ Ml. Cacciato et al. (2008) constrained the free 
parameters, criogz,, 71, 72, Mi and Lp, by fitting: the SDSS 
luminosity function of lBlanton et al.l ( 20031 ) and the galaxy- 
galaxy corre l ation lengths as a function of luminosity from 
IWang et al.l ([20071). The resulting best-fit parameters are 
listed in the first row of Table 1, and constitute our fidu- 



cial model Gl. We also consider two alternative models for 
$c(i|A/), called G2 and G3, the parameters of which are 
also listed in Table 1. For $s{L\M) we adopt the model of 
Cacciato et al. (2008) throughout, without any modifica- 
tions: i.e. models Gl, G2, and G3 only differ in P{AI\Lc) 
and have the same na^t{r\M). 

Having specified all necessary ingredients, we now com- 
pute the satellite weighted and host-weighted satellite kine- 
matics for our fiducial model Gl using Eqs. (|20|) and (|2H) . 
The results are shown as solid and dotted lines in Fig. [TJ 
where we have adopted a minimum satellite luminosity of 
imin = 10^ h~^ Lq. At the faint-end, the velocity disper- 
sions (Tsw and (Thw are equal, this simply reflects the fact 
that (A'aat)jvf —> ii Lc —> I/min- At the bright end, though, 
the non-zero scatter in P{M\Lc) causes the difference be- 
tween (Tsw and (Thw to increase systematically with increasing 
Lc- This is a generic trend for any realistic halo occupation 
model (see also van den Bosch et al. 2004). 

The upper panels of Fig. [5] show the mean and scatter 
of the MLR of central galaxies in models Gl (solid lines), 
G2 (dotted lines) and G3 (dashed lines). Panel (a) plots 
{logLc){M) — log(L*), which reveals the double power- law 
behavior of Eq. (|28|) . panel (b) shows the inverse relation, 



{logA/)(Le 



P{M\Lc)logMdM. 



(29) 



and panel (c) shows the scatter in the MLR, (TiogA/(Lc), 
deduced by using 
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0"lo 



P(A/|Lc)[logAf-(logA/)(ic 



dM . 



(30) 



Note that aiogM(Lc) increases with increasing Lc, even 
though the scatter aiogL is constant with halo mass. This 
simply owes to the fact that the slope of {logLc}{M) be- 
comes shallower with increasing Lc, as illustrated in Fig. |3l 

The comparison between models Gl and G2 illustrates 
the effect of changing the scatter aiogL in <E>c(i|A/). Both 
models have exactly the same (logLc)(Af) (the solid line 
overlaps the dotted line in panel a). However, because the 
scatter aiogL in G2 is larger than in Gl (see Table 1), the 
(log Af)(Lc) of G2 is significantly lower than that of Gl at 
the bright end (~ 0.5 dex at the bright end). This is due to 
the shape of the halo mass function. Increasing the scatter 
adds both low mass and high mass haloes to the distribution 
P{M\Lc) (cf. Eqs. [25] and [29]), and the overaU change in 
the average halo mass depends on the slope of the halo mass 
function. Brighter galaxies live on average in more massive 
haloes where the halo mass function is steeper. In particular, 
when the halo mass range sampled by P{M\Lc) lies in the 
exponential tail of the halo mass function, an increase in 
the scatter adds many more low mass haloes than massive 
haloes, causing a shift in the average halo mass towards 
lower values. On the other hand, fainter galaxies live in less 
massive haloes, where the slope of the halo mass function 
is much shallower. Consequently, a change in the scatter 
does not cause an appreciable change in the average mass. 
Finally, as expected, the scatter in the MLR, fJiogM{Lc), in 
G2 is higher than for Gl at all luminosities (see panel c). 

Panels (d) and (e) of Fig. [2] show the analytical predic- 
tions for (Jsw(ic) and (Jhw(ic), respectively. Note that mod- 
els Gl and G2 predict satellite kinematics that are signifi- 
cantly different (which can be distinguished given the typical 
measurement errors in Paper II) , even though both have ex- 
actly the same (logZ/c)(Af). In particular, model G2 predicts 
larger aaw and Uhw at the faint end, but lower osv and ahw at 
the bright end. The trend at the faint is due to the fact that 
the scatter o-iogj\/(ic) is higher in G2 than in Gl. Quantita- 
tively, this is evident from Eqs. (|10|l and (|lip . which demon- 
strate that both the satellite weighted and host weighted 
satellite kinematics increase with increasing scatter. At the 
bright end, however, the drastic decrease in {logM}{Lc) for 
G2 with respect to Gl overwhelms this boost and causes 
(Tsw and CThw to be lower in G2. 

Now consider model G3. This model has the same 
amount of scatter as model Gl, but we have tuned its param- 
eters (71, 72, Mi,Lq) that describe (log Lc}{M) such that its 
o"sw(ic) closely matches that of model G2 (the dotted and 
dashed curves in panel (d) are almost overlapping). As is 
evident from panels (a)-(c), though, the MLR of G3 is very 
different from that of G2. Note that the higher values of 
(logAf)(Lc) for G3 are compensated by its lower values of 
o'logAiiLc), such that the satellite- weighted kinematics are 
virtually identical. This clearly illustrates the degeneracy 
between the mean and the scatter of the MLR: One can de- 
crease the mean of the MLR and yet achieve the same ctsw 
by increasing the scatter of the MLR. It also shows that 
(Tsw alone does not yield sufficient information to uniquely 
constrain the MLR. 

Note, though, that although CTsw is the same for mod- 
els G2 and G3, their host-weighted satellite kinematics, 
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Figure 3. Illustration of the MLR of central galaxies. The solid 
black line indicates the mean of the Lc-M relation, while the gray 
scale region reflects the scatter. In this particular case the scatter 
in P{Lc\M) (indicated by vertical arrows) is taken to be constant 
with halo mass. Note, though, that the scatter in P(Af |Lc) (in- 
dicated by horizontal arrows) increases with increasing Lc; this 
simply is due to the fact that the slope of the mean Lc-M relation 
becomes shallower with increasing halo mass. 



o"hw(i'c), are different at the bright end. In fact, the ratios 
o"sw/o"hw for models G2 and G3 are clearly different. The 
logarithm of this ratio, shown in panel (f), follows the same 
trend as criog m (Lc) , i.e. it is higher for model G2 than for G3. 
This is in agreement with our toy model, according to which 
the ratio o-gw/fhw increases with the scatter criogAiiLc) (cf. 
Eq. [H]). This illustrates once again that the combination 
of (Taw and (Thw allows one to constrain both the mean and 
the scatter of the MLR simultaneously. 



6 SUMMARY 

The kinematics of satellite galaxies is a powerful probe of the 
masses of the dark matter haloes surrounding central galax- 
ies. With the advent of large, homogeneous redshift surveys, 
it has become possible to probe the mass-luminosity rela- 
tion (MLR) of central galaxies spanning a significant range 
in luminosities. Unfortunately, since most centrals only host 
a few satellite galaxies with luminosities above the flux limit 
of the redshift survey, one generally needs to stack a large 
number of central galaxies within a given luminosity bin and 
combine the velocity information of their satellites. Because 
of the finite bin-width, and because the MLR has intrinsic 
scatter, this stacking results in combining the kinematics of 
satellite gala:xies in haloes of different masses, which com- 
plicates the interpretation of the data. Unfortunately, most 
previous studies have ignored this issue, and made the over- 
simplified assumption that the scatter is negligible. 

Using realistic models for the halo occupation statis- 
tics, and taking account of selection effects, we have demon- 
strated a degeneracy between the mean and the scatter of 
the MLR; one can change the mean relation between halo 
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mass, M, and central galaxy luminosity, Lc, and simulta- 
neously change the scatter around that mean relation, such 
that the observed satellite kinematics, {(Tsat)(I/c), are unaf- 
fected. 

We have also presented a new technique to break this 
degeneracy, based on measuring the satellite kinematics us- 
ing two different weighting schemes: host-weighting (each 
central galaxy gets the same weight) and satellite weight- 
ing (each central galaxy gets a weight proportional to its 
number of satellites). In general, for central galaxies close 
to the magnitude limit of the survey, the average number of 
satellites per host is close to zero, and the satellite-weighted 
velocity dispersion, CTsw, is equal to the host-weighted veloc- 
ity dispersion, (Jhw This is because only those centrals with 
at least one satellite are used to measure the satellite kine- 
matics. For brighter centrals, however, (Taw > Thw and the 
actual ratio of these two values is larger for MLRs with more 
scatter (see Eq. [12] and panels c anf f of Fig. [2]). Hence, the 
combination of (Ts„(Z/c) and a"hw(ic) contains sufficient in- 
formation to constrain both the mean and the scatter of the 
MLR of central galaxies. In our companion paper (More et 
al. 2008), we apply this method to the SDSS, and show that 
the amount of scatter inferred from the data is in excellent 
agreement with other, independent constraints. In the com- 
panion paper, we also address the issues of measurement 
errors, sampling effects and interloper s. 

In a recent study, iBecker et al.l (|2007n analv zed the 
kinematics of MaxBCG clusters ( Koester et al.ll2007l ) and in- 
ferred the mean and the scatter of the mass-richness relation 
(here richness is a measure for the number of galaxies that 
reside in the cluster) . Becker et al. combined the kinematics 
of satellite galaxies in finite bins of cluster richness and mea- 
sured the second and fourth moments of the host-weighted 
velocity distribution. They used these two moments simul- 
taneously to determine the mean and the scatter of the 
mass-richness relation. This method is complementary to 
that presented here, and it will be interesting to compare 
both methods and investigate their relative strengths and 
weaknesses. We intend to address this in a future study. 

Finally we emphasize that the scatter in the condi- 
tional probability function P{M\Lc) is expected to increase 
with increasing Lc. This is due to the fact that the slope of 
(Lc)(M), which is the mean of P{Lc\M), becomes shallower 
with increasing halo mass. Hence, when stacking haloes ac- 
cording to the luminosity of the central galaxy, one can- 
not ignore the scatter in M, even when the scatter in 
P{Lc\M) is small. This has important implications for any 
technique that relies on stacking, such as satellite kine- 
matics and galaxy-g alaxy lensing (see e.g. iTasitsiomil 12004 : 
ICacciato et al.ll2008l) 
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